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Abstract 

The cooperative dynamics of a 1-D collection of Markov jump, inter¬ 
acting stochastic processes is studied via a mean-field approach. In 
the time-asymptotic regime, the resulting nonlinear master equation 
is analytically solved. The nonlinearity compensates jumps induced 
diffusive behavior giving rise to a soliton-like stationary probability 
density. The soliton velocity and its sharpness both intimately depend 
on the interaction strength. Below a critical threshold of the strength 
of interactions, the cooperative behavior cannot be sustained leading 
to the destruction of the soliton-like solution. The bifurcation point 
for this behavioral phase transition is explicitly calculated. 
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1 Introduction 

Interacting stochastic agents are modeled by a collection of nonlin¬ 
ear ly coupled Markovian stochastic processes. Inspired by the dy¬ 
namics recently exposed in [Balazs 2014], we focus on pure, right- 
oriented jump processes. For large and homogeneous swarms, the 
mean-field description offers a powerful method to characterize the 
resulting nonlinear global dynamics. Adopting the MF approach, the 
swarm behavior is summarized into a field density variable obeying 
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a nonlinear master equation. Such partial differential integral equa¬ 
tions are in general barely completely solvable. Nevertheless, several 
explicitly solvable models have been recently studied [Hongler 2014, 
Balazs 2014]. Our present goal is to enrich this yet available collec¬ 
tion by proposing an intrinsically nonlinear extension of the recent 
models introduced by [Balazs 2014]. Models involving pure jumps 
complete the solvable models with dynamics driven either by Brown¬ 
ian Motion and or by alternating Markov renewal processes [Hongler 
2014]. For strong enough mutual interactions, we explicitly observe 
the existence of a stationary probability measure propagating like a 
soliton. This soliton-like dynamics can be formed since the underly¬ 
ing nonlinear mechanism due to interactions exactly compensates the 
jump induced diffusion. This exhibits a close analogy with nonlin¬ 
ear wave dynamics where nonlinearity compensates the velocity dis¬ 
persion. Since the model is uni-dimensional, long-range interactions 
between the agents are mandatory for the existence of cooperative be¬ 
haviors here described by soliton-like probability measures. Decreasing 
the strength of the mutual interactions, via a barycentric modulation 
function similar to the one used in [Balazs 2014], we reach a critical 
threshold below which no stable cooperative behavior can be sustained. 
The critical threshold where the behavioral phase transition occurs can 
here be exactly calculated. 


2 Linear pure jump stochastic processes 

Let us first describe the dynamics of a single, isolated jump process 
which later in section 2, will enter into the composition of our inter¬ 
acting swarm. On M, we consider the right-oriented jump Markovian 
process X{t) characterized by the (linear) master equation: 

dtP{x,t) =-P{x,t) + f P{y,t)ip{x-y)dy, (1) 

J —OO 

where P{x,t) with P{x,0) = f{x) stands for the transition probability 
density. The function (p{x) : M —>• defines the probability density 

for the (right oriented) lengths of the process jumps. 

Taking the x-Laplace transform of Eq.([TJ and taking into account of 
the convolution structure, we obtain directly: 

dtP{s, t) = - [1 - (^(s)] P{s, t). (2) 

which solution reads: 

P{s,t) = (3) 
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where in writing Eq.([3|), we already did assume the initial condition: 


P{x,t) |t=o= 6{x). 


(4) 


Example. Consider the dynamics obtained when ^p{x) = \e~^^ yield¬ 
ing (p{s) = and when the initial probability density is f{x) = 6{x). 
Accordingly Eq.([3]) reads: 


P(s, t) = e ^ 



A+s > 
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The Laplace inversion of Eq.([5]) yields: 


(5) 


P{x, t) = e * < 


6{x) 




n—1 


^ re! (re — 1)! 
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:=J(x^t) 
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For J{x,t), we can write: 


J{x, t) 


A. 

dx 
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where Imiz) stands for the rre-modified Bessel’s functions. 
Hence the final probability density P{x,t) reads: 


(7) 


P{x, t) = e ^ 


6{x) + e 


—Xx 



X G M+ (8) 


and one may explicitly verify that one indeed has; P{x,t)dx = 1, 
(use the entry 6.643(2) in [Gradshteyn 80]). 

For time asymptotic regimes, Eq.® behaves as: 


lim P{x,t) ~ (9) 

t^QO 2y/Px't 

exhibiting therefore a diffusive propagating wave with vanishing am¬ 
plitude and velocity V ■= j- Due to translation invariance of the 
dynamics, we note that P(x — y,t) fulfills a 6{x — y) initial condition. 
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Hence, when P{x,0) = f{x), the linearity of the dynamics Eq.([T|) en¬ 
ables to write: 


Pf{x,0) = f{x), 

Pf{x,t) =4+ P{{x -y),t)f{y)dy. 


( 10 ) 


3 Non-linear Markovian jump processes 

Keeping the jumps probability density as (p{x) = let us now 

consider a large homogeneous collection of identical processes evolving 
like Eq.([T| now subject of mutual long-range interactions. The class of 
interactions we consider yields, in the mean-field limit, the nonlinear 
master equation: 


Q.{x,t) = f^g(z - (X(t)}) dzG{z,t)dz 
< dxtG{x, t) = -Q{x, t)da;G{x, t) + fl(y, t)dyG{y, t)Xe~P^~yUy, 

. = lM+ydyG{y,t)dy, 

( 11 ) 

where G{x, t) stands for the cumulative distribution of the a nonlinear 
jump process, (i.e. G{x, t) is monotonously increasing with boundary 
conditions G{—oo,t) = 0 and G{oo,t) = 1). Note that while in Eq.([T]) 
the jumping rate is unity, in Ea. (|lip it is replaced by Q{x,t) > 0 
which is explicitly state-dependent. This is precisely where the mutual 
interaction introduce a strong nonlinearity into the dynamics. In the 
sequel, we focus on cases where g{x) = g{—x) > 0. 

For asymptotic time, we now postulate that Ea. (|lip admits ^-functional 
dependent solutions with ^ = [x — Vt) and with the even symmetry: 

[ ^d^Gim = 0 , ( 12 ) 

JR 

where V is a propagating velocity parameter. In terms of Ea. dllll 
can be rewritten as: 

V [54G(0 + A4G(0] = 5^ . (13) 

Defining T(^) := log [cIjG(^)], after one integration step where the 
integration constant is taken to be zero, Eg. (1131) can be rewritten as: 

poo 

Vdl^CiO = -XV+ [g{g)dr,G{g)] dg. (14) 
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Assuming now a functional dependence g{^) = cosh“”(,^) with n € R, 
by direct substitution, it is immediate to see that Ea. (ll4l) is solved by 
the (normalized) probability density d^G{^): 


/ -pf m+1 \ 

> 0 , 

< m = X = 2 — n, (15) 


V = 


r(^) 

V^r(f)- 


Due to the ^-symmetry of the probability density Ea. (ll2l) is 

trivially fulfilled. 

For n £]2, —oo], Eo. dlbp implies that a stationary propagating density 
dxG{x) is sustained by the nonlinear dynamics Ea. (|lip . However, for 
short decaying 5 r(x)-modulation, occurring when n > 2, no stationary 
propagating probability density exists, (i.e. for this parameter range, 
m < 0 in Ea. (|15p and the solution cannot be normalized to unity as re¬ 
quired for a probability measure). For this exactly solvable dynamics, 
we also observe that the average jump length and the barycentric 
modulation strength controlled by the factor n are intimately depen¬ 
dent control parameters. In addition, we note that for large m, the 
asymptotic expansion of the F-function implies that lim^^cKD V — 


Illustration. Along the same lines as exposed in [Hongler 2014], the 
nonlinear dynamics given by Eq. m can be viewed as representing the 
mean-field evolution associated with a large population of stochastic 
jumping agents subject to a mutual imitation process. The swarm 
dynamics is described via the probability density function dxG{x, t) 
obeying a nonlinear partial differential equations (PDE). Agents mu¬ 
tual interactions are responsible for the state-dependent jumping rate 
D(x,t) in Eg. (HID . The functional form of Xl{x,t) simultaneously en- 
globes two distinct nonlinear features, namely: 


a) imitation process. To isolate this process, we may consider 
the case g{x) = 1, (i.e. n = 0) implying that 

Xl{x,t) = 1 — G{x,t). (16) 

The resulting state-dependent jumping rate Ea. dlbl) induces a 
traveling and compacting tendency. As the agents are subject to 
pure right-oriented jump, Ea. (ll6D effectively describes situations 
where the laggard agents jump more frequently than the leaders, 
(i.e. laggards try to effectively imitate the leaders behavior). 
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b) barycentric range modulation of the mutual interac¬ 
tions. The modulation obtained when g{x) 7 ^ 1 describes the 
relative importance attributed to interactions with agents remote 
from the barycenter {X{t)) of the swarm. Here, we may separate 
two distinct tendencies: 


i) when n G [ 0 , 2 [, far remote agents tend to not influence 
the dynamics. In this case, the resulting behavior can be re¬ 
ferred as a weak cooperate identity and the propagating 
probability density given by Ea. f|15D exhibits the shape of a 
table-top soliton with a plateau increasing when the lim¬ 
iting value 2 is approached. AOn observes a comparatively 
low propagating velocity V of these table-top like aggregates. 
Again, we emphasize that for n > 2, the cooperative interac¬ 
tions are not strong enough to sustain the propagation of a 
cooperative behavior in asymptotic time. This is well known 
in general for 1-D stochastic interacting system, (the Ising 
model being the paradigmatic example) where no coopera¬ 
tive phase can be formed when the interactions operate on 
too limited ranges. 


ii) for n < 0 , the g(x) modulation effectively gives rise 
to a strong cooperate identity. Far remote agents in¬ 
creasingly influence the swarm. This gives rise to sharply 
peaked solitons-like probability densities propagating with 
high propagating velocities. 
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